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Abstract. We characterize the crossover regime to the KPZ equation for a 
class of one-dimensional weakly asymmetric exclusion processes. The crossover 
depends on the strength asymmetry an^~^ (a, 7 > 0) and it occurs at 7 = 1/2. 
We show that the density field is a solution of an Ornstein-Uhlenbeck equation 
if 7 G (1/2, 1], while for 7 = 1/2 it is an energy solution of the KPZ equation. 
The corresponding crossover for the current of particles is readily obtained. 



1. Introduction 

One-dimensional weakly asymmetric exclusion processes arise as simple 
models for the growing of random interfaces. For those processes the micro- 
scopic dynamics is given by stochastic lattice gases with hard core exclusion 
with a weak asymmetry to the right. The presence of a weak asymmetry 
breaks down the detailed balance condition, which forces the system to exhibit 
a non trivial behavior even in the stationary situation. Using renormaliza- 
tion group techniques, the dynamical scaling exponent has been established 
as z = 3/2 and one of the challenging problems is to derive the limit distribu- 
tion of the density and the current of particles Il23l . 

For asymmetric exclusion processes, partial answers have been given in 
particular settings as starting the system from the stationary state and from 
specific initial conditions. For these models, under a certain spatial shifting 
and time speeding, the current of particles has Tracy- Widom distribution, see 
IfTTl [I6l [24l [25 1 ■ The Tracy- Widom distribution was initially obtained in the 
context of large N statistics of the largest eigenvalue of random matrices, but 
has been recently obtained as the scaling limit of stochastic fields of random 
models, see HI [6l [I6l HH |22 1 and references therein. 

Here we are interested in establishing the equilibrium density fluctuations 
for weakly asymmetric exclusion processes with strength asymmetry an?^'' , 
that we fully describe below We consider the system under the invariant state: 
a Bernoulli product measure of parameter p e [0, 1] that we denote by Vp. By re- 
lating the current with the density of particles, as a consequence of last result 
we derive the equilibrium fluctuations of the current of particles. 

The weakly asymmetric simple exclusion process was studied in (Tj [9j, for 
strength asymmetry n (that corresponds to 7 1 in our case), and in ['51 for the 
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strength asymmetry r?!-'' (that corresponds to 7 = 1/2 in our case). For 7 = 1, 
the equihbrium density fluctuations are given by an Ornstein-Uhlenbeck pro- 
cess which imphes the current to have Gaussian distributions, see (ZllQj. For 
7 = 1/2, fSl used the Cole-Hopf transformation to derive the non-equihbrium 
fluctuations of the current. By using the Cole-Hopf transformation initially 
introduced in Il8), one obtains an exponential process. For this exponential pro- 
cess, the limiting fluctuations are given by the stochastic heat equation, which 
is a linear equation, making the asymptotic analysis much easier. As a con- 
sequence, the fluctuations of the original process are easily recovered. In our 
approach, we study the weakly asymmetric exclusion directly, but in order to 
identify the limiting density field as a weak solution of a stochastic partial dif- 
ferential equation, we have to overcome the difficulty of closing the equation 
by means of the Boltzmann-Gibbs principle. For this reason our results are 
restricted to the equilibrium setting. 

It is known from |7, 9| that for 7 = 1, the limit density field is a solution of 
an Ornstein-Uhlenbeck equation which has a drift term. This drift term comes 
from the asymmetric part of the dynamics and can be removed by taking the 
process moving in a reference frame with constant velocity. By removing the 
drift of the system, there is no effect of the strength of the asymmetry on the 
distribution of the limit density field. In order to see how far this picture can 
go, we strengthen the asymmetry by decreasing the value of 7. We can show 
that, for 7 e (1/2, 1] there is still no effect of the strength of the asymmetry on 
the limiting density field. In this case, the limiting density field is still solution 
of the Ornstein-Uhlenbeck equation as for 7 = 1 and for this reason the process 
belongs to the Edwards- Wilkinson IflOil universality class. Nevertheless, for 
7 = 1/2 the limiting distribution "feels" the effect of the strengthening of the 
asymmetry, by developing a non linear term in the equation that characterizes 
the limiting density field. In this case the limiting density field is a solution 
of the Kardar-Parisi-Zhang (KPZ) equation, so that for 7 = 1/2 the process 
belongs to the KPZ [17 1 universality class. 

The KPZ equation was proposed in iflTll to model the growth of random in- 
terfaces. Denoting by ht the height of the interface, this equation reads as 

dth = DAh + a{Vhf + aWt, 

where D, a, a are related to the thermodynamical properties of the interface 
and Wt is a Gaussian space-time white noise with covariance given by 

E\Wt{u)-WsW)\ = 8{t - s)6(u - v). 

According to the dynamical scaling exponent z = 3/2, a non-trivial behavior 
occurs under the scaling hn{t,x) = n^^/^h{tn'^/^ ,x/n). This means, roughly 
speaking, that in our case, for 7 = 1/2 a non trivial behavior is expected even in 
the stationary situation and in that case the model belongs to the universality 
class of the KPZ equation. 

To our knowledge, a rigorous mathematical proof of the characterization 
of the intermediate state between the Ornstein-Uhlenbeck process and the 
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crossover to the KPZ equation, was lacking so far, see IfTSlI and references 
therein. Nevertheless we refer the reader to the paper IH in which the authors 
characterize the crossover regime for a special case of weakly asymmetric ex- 
clusion process, but going through the Cole-Hopf transformation. 

As a consequence of obtaining the fluctuations of the density of particles, 
we obtain the equilibrium fluctuations for the current of particles for different 
strength regimes depending on the strength asymmetry. More precisely, we 
show that for 7 > 1/2 the current properly centered and re-scaled converges to 
a fractional Brownian motion with Hurst parameter 1/4 and for 7 = 1/2 the 
limit process is given in terms of the solution of the KPZ equation. 

The existence of a non-trivial crossover regime for weakly asymmetric sys- 
tems was found in |2|. There, the authors use the theory developed in |3| and 
m, and found that for a phase of weak asymmetry the fluctuations of the cur- 
rent are Gaussian, while in the presence of stronger asymmetry they become 
non-Gaussian. Here, we provide the characterization of the transition from 
the Edwards-Wilkinson class to the KPZ class, for general weakly asymmetric 
exclusion processes. We prove that the transition depends on the strength of 
the asymmetry without having any other intermediate state and by establish- 
ing precisely the strength in order to have the crossover. We point out here, 
that our results are also valid for weakly asymmetric exclusion processes with 
finite-range interactions. All the proofs follow in this case with minor nota- 
tional modifications. 

Here follows an outline of this paper. In the second section, we introduce 
the model and we describe the equilibrium density and current fluctuations for 
the process under different strength asymmetry regimes. In the third section, 
we sketch the proof of the results for the intermediate state regime and in 
the fourth section we recall briefly the results about the crossover to the KPZ 
equation established in Ifl3]| . 

2. Equilibrium fluctuations 

Let T]t be the weakly asymmetric exclusion process evolving on Z. The state 
space of this Markov process is Q {0, 1}^ and its dynamics can be described 
as follows. On a configuration 77 e and after a mean one exponential time, a 
particle jumps to an empty neighboring site according to a transition rate that 
has a weak asymmetry to the right and depends on a function c(ry). We assume 
c : f2 — >^ M to be a local function, bounded from above and below, and that turns 
the system gradient and reversible with respect to the stationary state Vp. The 
gradient condition is the most restrictive one and requires the existence of a 
local function : — > M such that for any 77 e 

The process is speeded up on the diffusive time scale so that 77" = r]tn2 . 
Here and in the sequel, for r; e and a; e Z we denote by t^ij the space trans- 
lation by X, namely for y G Z, Tx7]{y) — r]{y + x) and for a function / : il — > M 
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we denote by Txf{ri) the induced translation in /, namely Txf{ri) ~ /(tx??). Let 
Cx{ri) = Txc{rj). For the configuration ?/, the transition rate from x to x + 1 is 
given by Cx{r])pn and fromx+l to x is given by Cx{r])qn, where p„ := (l+a/n'^)/2, 
Qn '■— I — Pn and a > 0. We refer the reader to flSl for a complete discussion 
about the assumptions on c(-). 

The generator of this process acts over local functions / : O ^ M as 

Ai/(?7) ^n'^^ c^x{v){P7iV{x){l - v{x + 1)) + QnVix + 1)(1 - 'n{x))}W^^x+if{r]), 

where n e N, Vx^x+if{r]) ^ fiv'''''^^) - fiv) and 

{r]{x + 1), z = X 
T]{x), z^x + 1 
ri{z), z ^ x,x + 1. 

If a = the process 77" is said to be symmetric and if c( ) = 1 the process is 
the symmetric simple exclusion process. On the other hand if c( ) = 1, a = 1 
and 7 = 1, the process is the weakly asymmetric simple exclusion process, 
studied in [T, 91, and that can be interpreted as having the symmetric and 
asymmetric dynamics speeded up by and n, respectively. 

We notice that decreasing the value of 7 in p„ above, corresponds to speeding 
up the asymmetric part of the dynamics on longer time scales as n'^^'*, while 
the scaling of the symmetric dynamics is not affected. 

A stationary state for this process is the Bernoulli product measure on O of 
parameter p e [0, 1] that we denote by i^p and whose marginal at ri{x) is given 

^y^piv ■ = 1) = p- 

2.1. Hydrodynamic Limit. Here we recall briefly the hydrodynamic limit 
for 77". The hydrodynamical scaling corresponds to the strength asymmetry n. 

For that purpose we introduce the empirical measure as the positive mea- 
sure in R defined by 

= - V77t"(a;)(5^/„(dx), 

xez 

where for u G M, (5„ is the Dirac measure at u. 

Take po : M — > [0, 1] a strictly positive and piecewise continuous function 
such that there exists p e (0, 1) satisfying / \p{){x) — p\dx < +00. Start the 
process 77" from {^„; 71 e N} - a product measure on fl, whose marginal at ri{x) 
is Bernoulli of parameter po{x/n), namely: 

^J■n{v ■ V{x) = 1) = Po{x/n). 

Then, TT^{dx) converges in probability to the deterministic measure p{t, x)dx, 
where {pit, x); i > 0, x e R} is the unique weak solution of the viscous Burgers 
equation 

dtp{t, x) = ^A^h{p{t, x)) - aV Pipit, x)) 
p(0, u) — poiu),u e M. 



CROSSOVER TO THE KPZ EQUATION 



5 



where iph{p) = J hdvp, l3{p) — x{p) J cdvp and xip) = pi^ - p)- We notice that, 
thanks to the gradient (reversible) condition of c( ), (ph{p) (respectively a/3(p)) 
corresponds to the expectation with respect to Vp of the instantaneous current 
of the symmetric (respectively asymmetric) part of the dynamics. Our purpose 
here is to analyze the fluctuations of the empirical measure from the stationary 
state Vp and from there, to derive the fluctuations of the current of particles. 

2.2. Equilibrium Fluctuations. From now on we fix a density p e (0, 1). De- 
note by 5(K) the Schwartz space, i.e. the space of rapidly decreasing functions 
and denote by S'{M) its dual with respect to the inner product of i^(R). Let 
2?([0,cx)),5'(R)) be the space oi cadVag trajectories from [0,oo) to 5'(M). De- 
fine {yt'it > 0} as the density fluctuation field, a linear functional acting on 
H e 5(M) as 

ynH)^^Y.H{l)m-)-p)- (2.1) 

We start by considering 7 = 1, which as mentioned above corresponds to the 
hydrodynamic strength asymmetry. By computing the characteristic function 
of [H) it follows that J^q converges in distribution to a spatial white noise of 
variance x{p) - the static compressibility of the system. Moreover, {3^"; n e N} 
converges in distribution with respect to the Skorohod topology of the space of 
cadlag trajectories I?([0, oo), 5'(M)) to the process yt, solution of the Ornstein- 
Uhlenbeck equation 



dyt = ^^Aytdt - ai3'{p)vytdt + .Mp)VdWt, (2.2) 

where Wt is a space-time white noise of unit variance. To prove this result, 
since the behavior at t = is characterized it remains to analyze the time 
evolution of the limit density field. For this purpose we introduce the martin- 
gale associated to 02. ID and we analyze its asymptotic behavior. The martin- 
gale decomposition presents two integral terms which cannot be written as a 
function of the density fluctuation field given in ( 12.11 1. As a consequence we 
cannot identify straightforwardly the limit density field as a weak solution to 
some stochastic partial differential equation. In order to perform this identi- 
fication, a replacement argument is required. This replacement is known as 
the Boltzmann-Gibbs principle and was introduced in Ifl9]| . With this result in 
hand and since the system is gradient, the equilibrium fluctuations are easily 
derived, see for example [71 for a detailed proof 

Now we analyze equation ( |2.2| l. This equation has a drift term that vanishes 
if P' (p) is equal to 0, and in that case the equation does not depend on a any 
longer. In order to remove the drift term (which arises from the asymmetric 
part of the dynamics) from the limiting density field, we take 77" moving in a 
reference frame with constant velocity al3'{p)n. In order to see the dependence 
on the strength of the asymmetry on the limit process, we redefine the density 
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fluctuation field on H e 5(M) as: 




(2.3) 



where T^H{■) = H{- — aj3' {p)tn^^^) and aP'{p)n^~^ is the velocity of the sys- 
tem (be aware of the scahng parameter x — > x/n). As above, it is not hard 
to show that for 7 = 1, {y^'^•,n S N} converges in the Skorohod topology of 
I?([0, oo),5'(]R)) to the process yt solution of the Ornstein-Uhlenbeck equation 



dyt = ^^^^ytdt + V^VdWt, (2.4) 

which corresponds to equation ( I2.2I > with a = 0. Since for 7 = 1 the fluctuations 
of the density are given by ( 12. 4D . the system belongs to the Edwards- Wilkinson 
universality class |10|. 

In order to see the effect of incrementing the strength asymmetry in the 
limit density field we decrease the value of 7. As discussed in ||5l, the effect of 
the asymmetry is presented in the limiting density field when 7 = 1/2 and in 
that case yt has a very different qualitatively behavior from the one obtained 
for 7 = 1, namely the solution of ( I2.4I >. Here we characterize the limiting 
density field yt for the intermediate state regime, namely we show that: 

Theorem 2.1. If ^ ^ (1/2,1], then the sequence {y^'^^-^n e N} converges in 
distribution with respect to Skorohod topology ofV^O, 00), 5'(M)) to the process 
yt solution of the Ornstein-Uhlenbeck equation ( I2.4I >. 

By the previous result, we have that for 7 e (1/2, 1] the system still belongs 
to the Edwards- Wilkinson universality class. So, if we want to see the effect 
of strengthening the asymmetry in the limiting density field, we have to take 
7 = 1/2 which is in agreement with the result of IS). Recently in Ifl3l . it was 
shown that for 7 = 1/2, the sequence {3^"'^; n g N} is tight and any limit point 
is an energy solution of the KPZ equation: 



dyt = ^Aytdt ~ !^li^v{y, fdt + Mp)^dWt. (2.5) 

So, our result says that weakly asymmetric simple exclusion processes be- 
long to the Edwards-Wilkinson universality class for 7 e (1/2,1]; and cross 
to the KPZ universality class when the strength asymmetry is precisely given 
by an3/2, see |13]| . 

3. Beyond the hydrodynamic time scale 

In this section we give an outline of the proof of the equilibrium density and 
current fluctuations in the intermediate state regime (with 7 e (1/2, 1]), i.e. 
between the Ornstein-Uhlenbeck process ( 12. 4D and the KPZ equation (|2]5). 
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3.1. Density fluctuations. Here we prove Theorem |27TJ Recall the definition 
of the density fluctuation field given in i2.3i . In this setting, we remove the 
drift of the process, so that we suppose to have r]" moving in a reference frame 
with constant velocity given by a(3' {p)n'^^''' with 7 e (1/2, 1]. 

At t — 0, by computing the characteristic function of 3^o this field con- 
verges to a spatial white noise of variance Now, we analyze the asymp- 
totic behavior of some martingales associated to 3^"'^, in order to identify the 
limit density field yt as a weak solution of the stochastic partial differential 
equation ( 12.41 1. For that purpose, fix iJ e S(R) and notice that by Dynkin's 
formula. 



is a martingale with respect to the natural filtration Ft = (j{ijs,s < t), where 

/(r;)=ac(,7)(,7(l)-77(0))V2, 

r^^Vfiv) - Tj{r^) - aPip) ~ aP'{p){r,{x) - p), 

A"iJ(a;/n) := n^{H{{x + l)/n) + H{{x - l)/n) - 2H{x/n)) and 
V"H{x) := n{H{{x + l)/n) - H{x/n)). 

Notice that A" and V" are the discrete Laplacian and the discrete derivative, 
respectively. We point out that the mean of / with respect to Vp is given by 

am- 

The quadratic variation of equals to 

We notice that if 7 > 0, the term corresponding to a/rf inside last integral, 
vanishes in L'^{¥^^) as n — > +00. From now on, Pi.^ denotes the distribution of 
the Markov process 77" starting from the stationary state Vp and Ej^^ denotes 
the expectation with respect to P^^ . Last term arises from the asymmetric part 
of the dynamics and for this reason, only for a big strength of the asymmetry 
(namely for 7 = 0; which corresponds to speeding up the asymmetric part of 
the dynamics by v?) it will give rise to additional stochastic fluctuations of the 
system. From simple computations, it follows that the limit as n +03 of 
the martingale M'^'^iH) is given by \\^/l3{p)V H\\2Wt{H), where Wt{H) is a 
Brownian motion and || • II2 denotes the L^(R)-norm. 
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Now, we need to analyze the limit of the integral terms. We start by the less 
demanding, namely Invoking the Boltzmann-Gibbs principle intro- 

duced in lfT9l . can be written as 



plus an i^(Pyp) negligible term. 

Now we analyze We notice that if we were considering the case 

7 = 1 the Boltzmann-Gibbs principle as stated in |19|, would be saying that 
A"''^{H) vanishes in (P^p ) as n +cx). Since7 < 1, the result in [19 1 does not 
give us the information we need about the limit of this integral term. Never- 
theless, according to Corollary 7.4 of [12|, it follows that in fact for 7 e (1/2, 1] 
the same result is true. Indeed the integral term Ai''^{H) still vanishes in 
i^(Pi,p) as n — > +00 for 7 e (1/2, 1]. We remark that the mentioned result in 
Ifl2]| was proved for the symmetric simple exclusion process but it is also true 
for our model. More explicitly, for our case that result says the following: 

Proposition 3.1 (Stronger Boltzmann-Gibbs Principle [121). 

Let tp : V, ^ M. be a local function and let <f^,{p) := Ei,^ [V-'('7)]- Por 7 e (1/2, 1] 
and H e S{W), it holds that 



lim 

n— ^00 



t 1—^ 

n ^ 



/n 



= 0. 



In fill, this result was proved for ip{ri) = (77(0) — p){ti{1) — p), but the proof 
holds for any local function, since the fundamental ingredients invoked along 
the proof, namely the spectral gap bound and the equivalence of ensembles, 
hold in general for local functions. As a consequence of last result, ^"'^(i/) 
still vanishes in L^(P^p) as n +00 for 7 e (1/2, 1]. 

Putting together the previous observations, for 76 (1/2, 1], the limit density 
field 3^t(i?) satisfies: 

yt{H) = yo{H) + J^ys{^^AH)ds + \\^^WH\\2Wt{H), 

so that yt is a weak solution of ( 12. 4D . With this decomposition it follows that 
the covariance of the limit field is given on H,G £ S{R) by 

E[yt{H)y,{G)] = xip) I Tt,sH{x)G{x)dx, (3.1) 
Jr 

where {ri}t>o is the semigroup associated to the operator ^^^^A. 

In order to complete the argument, it remains to show tightness of the se- 
quence {y^'^^jneN- In lfT2ll this was shown for the asymmetric simple exclusion 
process, but the same computations with minor modifications, hold for the pro- 
cesses we consider here. 
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3.2. Current fluctuations. Now we want to derive the fluctuations of the 
current of particles, from the fluctuations of the density of particles. Since we 
took the process moving in a reference frame with constant velocity a(3'{p)n'^^'^ , 
we consider the current of particles through a moving bond. 

Fix 7 e (1/2, 1]. For a site x denote by J^it), the current of particles through 
the bond {x, x + 1}, i.e. the number of particles that jump from the site x to 
x + 1, minus the number of particles that jump from x + ltox during the time 
interval [0,tn^]. Formally 

y>x~\-l 

so that it is the difference between the density fluctuation field ( 12. Il l at time t 
and at time 0, evaluated on the Heaviside function = 1(2;, oo)- 

Consider the line — x + [a/3' {p)tn^^'^] and let J^''^{t) be the current of 
particles through the time-dependent bond {ax-, a^ + 1}. For a; e M, [x] denotes 
the integer part of x. Take a; = to simplify the exposition but for any other 
site the results stated below are also true. 

Up to the strength asymmetry n'^~^ with 7 > 1/2, we are able to show that 
the current properly centered and re-scaled converges to a fractional Brownian 
motion of Hurst parameter H = 1/4: 

Theorem 3.2. Fix a; e Z, 7 e (1/2, 1] and let 

Then, for every k > 1 and every < ti < t2 < ■■■ < t^, (Zj" , ZJ^J converges in 
law to a Gaussian vector {Zt^ , ) with mean zero and covariance given by 

E[ZtZs] = \l^^^x{p){Vt + - Vt^) 

provided s <t. 

We sketch here the proof of last result. The idea of the argument is to obtain 
the fluctuations of the current from the fluctuations of the density of particles, 
namely from Theorem l2.1l This argument was initially proposed in ||20| . For 
more details we refer the reader to 112] and II15I . 

From the stronger Boltzmann-Gibbs principle as stated in Proposition 13.11 
it follows that for 7 e (1/2,1] 

vanishes in L^{F^^) as £ ^ +(X) uniformly over n, where {Gejeen is a sequence 
approximating the Heaviside function Hq, which can be taken for example 
equal to Gi{x) = (1 - x/£)+. 

Combining last result with the convergence of 3^"'^, it follows that con- 
verges to a random variable, which formally reads as 
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where yt is the solution of the Ornstein-Uhlenbeck equation ( I2.4I I. 

The same argument can be apphed to show the same result for any vector 
(Ztj , Zt^). To compute the covariance we do the following: 

E[ZtZ,] = E[{yt{Ho) - yoiHo)}{ys{Ho) - yoiHo)}] 

= lim E[{ytiGi)~yo{Ge)}{ys{Ge)~yo{Ge)}]. 

l—y+oo 

Now we use 03.11 1 to write last expression as 

X{p) lim / (Tt-,Gi{x)Gfix) - Ge{x)TtGiix) - GeXx)T,Gi{x) + G^,{x))dx. 

Using the definition of Tt{Gi) we get the covariance stated in the theorem. 

We notice that last convergence takes place in the sense of finite-dimensional 
distributions. Since the distributions of >'t(i7o) are Gaussian, this implies the 
limit current to be Gaussian distributed. 

4. The Crossover regime 

Here we describe briefly the limit density fluctuation field for 7 = 1/2, by 
recalling the arguments used in [.13.1 . As a consequence, the fluctuations of the 
current are easily derived. 

4.1. Density fluctuations. Recall the definition of the density field given in 
( 12.31 1 with 7 = 1/2. We want to obtain the limiting density field as a solution of 
the stochastic partial differential equation ( 12.51 1. 

For that purpose we introduce the martingales associated to 3^"'^ as in the 
previous section. The main difference on the limiting density field for this 
strength of the asymmetry, comes from the limit of the integral term At'''{H). 
We saw above that if 7 > 1/2 the somehow "stronger" Boltzmann-Gibbs princi- 
ple derived in |12|, tells us that A^'^'iH) vanishes in L^{P^^) as n +00. 

In the presence of a stronger asymmetry, a second order Boltzmann-Gibbs 
principle is needed and in [13J it was derived. This stronger replacement is 
derived through a multi-scale argument that was introduced in II 121 . This 
multi-scale approach, allows to obtain the desired replacement as a sequence 
of minor replacements in microscopic boxes that duplicate size at each step, 
until a point in which the sum of the errors committed at each step is neg- 
ligible and the last replacement holds at a box of small macroscopic size. In 
this macroscopic box, the non-trivial term can be identified as the square of 
3^"'^- For details on this argument we refer the reader to [ 13] . The funda- 
mental features of the model that are used in order to derive this second order 
Boltzmann-Gibbs principle, are the sharp spectral gap bound for the dynamics 
restricted to finite boxes, plus a second order expansion on the equivalence of 
ensembles. These results are quite general and this argument can be applied 
for more general models than of exclusion type. 
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With this procedure, in |[T3]| was shown that for any H e S{R), A'l''{H) 
converges in {P^p ) to 

MH) = lim-^^ /* / y,{i,{x)fH'{x)dxds. (4.1) 

Here ie{x){y) = e^^l(x<y<x+^)- 

Applying the same arguments as above, in IfTSl it was shown that any hmit 
point of 3^"'^ satisfies: 

yt{H) - y^iH) - J^y,(^^^AH)ds + AtiH) + \\y/Mp)'^H\\2WtiH) 

with AtiH) given as in ( I4.1I ). According to [13|, the Umit density field yt is a 
weak solution of the KPZ equation ( 12. 5D . 

4.2. Current fluctuations. As mentioned above, having established the fluc- 
tuations of the density of particles one can obtain the fluctuations of the cur- 
rent of particles J^^'^{t). Following the route described above, it follows that 

-l={j';'^(i)-E,^[j';'^(i)]} 

converges to 

yt{H^) - yo{Ho) 

where yt is the limit of the subsequence y^ and it is a weak solution of the 
KPZ equation ( l23] l. 

So, the crossover regime for the current of particles occurs from Gaussian 
to a distribution which is given in terms of the solution of the KPZ equation 
(|2]5]l. As for the density of particles, the crossover of the current also occurs at 
strength asymmetry ar?!'^ . 

In [1, 22 1 the authors studied the weakly asymmetric simple exclusion pro- 
cess (which corresponds to taking c( ) = 1 here) starting from different initial 
conditions. There, the crossover regime for the current of particles is studied 
and the transition goes from Gaussian to Tracy- Widom distribution. For the 
strength asymmetry r?^"^ , |22 | identify the limit distribution of the current as 
the difference of two Fredholm determinants which, from the random matrix 
theory, is known to converge to the Tracy- Widom distribution. 

In [131, the stationary situation is studied and the limit of the current is 
written in terms of the energy solution of the KPZ equation. The result is 
true for a general class of weakly asymmetric exclusion processes. Here we 
show that the transition from the Edwards-Wilkinson class to the KPZ class is 
universal within this class of processes and the strength asymmetry in order 
to have the crossover is precisely ar?!"^ . This is a step towards characterizing 
the universality of the KPZ class. 
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